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q \ Abstract: The present work shows that through a suitable change of variables rela- 

tivistic dynamics can be mapped to light propagation in a non-homogeneous medium. 
A particle's trajectory through the modified space-time is thus formally equivalent to a 
light ray and can be derived from a mechanical equivalent of Fermat's principle. The 
similarities between light propagation and mechanics are then extended to quantum me- 
chanics, showing that relativistic quantum mechanics can be derived from a wave equa- 
tion in modified space-time. Non-relativistic results, such as de Broglie's wavelength, 
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Schrodinger equation and uncertainty principle are shown to be direct consequences of 
the theory and it is argued that relativistic conclusions are also possible. 

1 Introduction 

This paper is presented in a rather crude state; the text is imperfect and some conclusions are 
deferred to ulterior publications; nevertheless the author feels that this work must be diffused even 
in a preliminary stage due to its significance. The author's presence at the OSA's annual meeting 
O ■ provided an opportunity for the presentation of his work that he could not despise. 

The similarities between light propagation and wave mechanics have been pointed out by numer- 
ous authors, although a perfect mapping from one system to the other has never been achieved. 
q ; Almeida et al. 1 showed that near- field light diffraction could be calculated using the Wigner Dis- 
tribution Function (WDF) and obtained results proving the existence of super-resolution in certain 
circumstances. 

. 

The study of wide angle light propagation makes use of a transformation which brings to mind the 
Lorentz transformation of special relativity. It was then natural to try an association of Newtonian 
^ ■ mechanics to paraxial optics and special relativity to wide angle propagation. This process pro- 
moted the definition of a coordinate transformation to render the relativistic space homologous to 
the optical space. The introduction of a variational principle allowed the derivation of relativistic 
dynamics in the modified space-time in a process similar to the derivation of optical propagation 
from Fermat's principle. One important consequence is that each particle travels through modified 
space-time with the speed of light. 

The similarity could be carried further to diffraction phenomena and quantum mechanics. It was 
postulated that a particle has an intrinsic frequency related to its mass and many important results 
were derived directly from this statement. More general results will probably be feasible in the 
future. 



Almeida, Optical interpretation. 



OSA Annual meeting/2008 Page 2 



2 Notes on Hamiltonian optics 
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The propagation of an optical ray is governed by Fermat's principle, which can be stated: 

5 [ P \ds = . (1) 
JPi 

The integral quantity is called point characteristic and measures the optical path length between 
points Pi and P 2 - 

V(x ll y 1 ,x 2l y 2 ) = nds . (2) 
J Pi 

The quantity ds is the length measured along a ray path and can be replaced by: 

ds = — = — (3) 

(i — u\ — u*yi 2 u z 

where u x , u y and u z are the ray direction cosines with respect to the x, y and z axes. 
Inserting Eq. (3) into Eq. (1) we have: 

rz 2 ndz 
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Of course we can also write: 
with 



ds = (1 + x 2 + y 2 ) 1/2 dz , (5) 



dx dy 
x = — and y — — . (6) 



dz dz 



It is easy to relate x and y to u x and u 
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where only the positive root is considered. 
Inserting Eq. (5) into Eq. (1) we get: 



5 f 2 n(l + x 2 + y 2 ) l/2 dz = . 



(8) 



We can use the position coordinates x and y as generalized coordinates and z for time, in order to 
define the Lagrangian. We have: 3 ' 4 

L = n(l + x 2 + y 2 ) 1 ' 2 (9) 
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Euler Lagrange's propagation equations are: 

d dL dL 

dz ox ox 

± d ±- d ± = 0; (11) 
dz dy dy 

We can go a step further if we define a system Hamiltonian and write the canonical equations; we 
start by finding the components of the conjugate momentum (p) from the Lagrangian. Knowing 
that n is a function of x, y and z, the conjugate momentum components can be written as: 

nx 



dL 

d± ~ {l + x 2 + y 2 y/ 2 
dL ny 
dy~ = (l + x 2 + y 2 y/ 2 ' 



Px = 

dL _ ny 

Py ~ a/. fl , A.2 , „-.2M/2 ' > 



If we consider Eq. (7), the result is: 



The system Hamiltonian is: 



Px TlUx j 

Py = UUy . (13) 



H = p x x+p y y-L 
n(x 2 + y 2 ) 



x 2 + y 2 ) 1 / 2 
—n 



(l + x 2 + y 2 y/ 2 
-n{l-ul-ulf/\ 

-nu z , (14) 



The Hamiltonian has the interesting property of having the dependence on the generalized co- 
ordinates and time, separated from the dependence on the conjugate momentum. The canonical 
equations are: 



dH u T u 



x 



x 



dp x (1 - u 2 x - uf) 1 / 2 u z 

8H Uy Uy 

dp y (1 - u 2 x - u 2 ) 1 / 2 u z ' 

dH 2 2 U7 dn dn 

nu x + nu x = -— = { l-u x -u y )l - = u z - , 

dH 2 2xi /9^ n 

nu y + nu y = -— = {l- Ux -u y ) — = u z — . (15) 

Obviously the first two canonical equations represent just a trigonometric relationship. 
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It is interesting to note that if the refractive index varies only with z, then the conjugate momentum 
will stay unaltered; the direction cosines will vary accordingly to keep constant the products nu x 
and nu y . 

We will now consider an non-homogeneous medium with a direction dependent refractive index 
and will add this dependence as a correction to a nominal index. 

n c 

n = n o~ . 1/2 , (16) 

(1 + x l + y l ) ' 

where n is the nominal index and n c is a correction parameter. Eq. (9) becomes 

L = n (l + x 2 + y 2 ) - n c . (17) 

We will follow the procedure for establishing the canonical equations in this new situation. It is 
clear that the momentum is still given by Eq. (13) if n is replaced by n . 

The new Hamiltonian is given by 

H = -n (l - u 2 x - u 2 ) + n c , (18) 



and the canonical equations become 



dH u r u 
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y 



n u x + n u x = — — = u z 



dp x (1 -u 2 x -u 2 y ) x l 2 u z 

8H Uy Uy 

dp y (i - u 2 x - u 2 y y/ 2 

dH dn n dn r 



X 
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dx dx dx 

. dH dn dn c 

n u y + n u y = — — = u z — — . (19) 

ay ay ay 

The present discussion of non-homogeneous media is not completely general but is adequate for 
highlighting similarities with special relativity and quantum mechanics, as is the purpose of this 
work. 



3 Diffraction and Wigner distribution function 

Almeida et al. 1 have shown that the high spatial frequencies in the diffracted spectrum cannot 
be propagated and this can even, in some cases, lead to a diffraction limit much lower than the 
wavelength; here we detail those arguments. 

The Wigner distribution function (WDF) of a scalar, time harmonic, and coherent field distribution 
<£>(q, z) can be defined at a z = const, plane in terms of either the field distribution or its Fourier 
transform <^(p) = / </>(q) exp(— z£;q T p)<iq: 5 ~ 7 

W(q,p) = /^(q+yj^^q-yjexp^q^dq' (20) 
= ^ J Tp (p + ^) p (p - ^) ex P (zfcq T p') dp' , (21) 
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where k = 2n/\, * indicates complex conjugate and 

q = (x,y), (22) 

p = (nu x ,nu y ). (23) 

In the paraxial approximation, propagation in a homogeneous medium of refractive index n trans- 
forms the WDF according to the relation 

W(q,p,z) =W(q--p,p,0). (24) 

n 

After the WDF has been propagated over a distance, the field distribution can be recovered by 6 ' 7 

cp(q, z)cp*(0, z) = J W(q/2, p, z) exp(zqp)dp. (25) 

The field intensity distribution can also be found by 

4tt 2 r 

\cp(d,z)\ 2 = — J W(q,p,z)dp. (26) 

Eqs. (25) and (26) are all that is needed for the evaluation of Fresnel diffraction fields. Consider 
the diffraction pattern for a rectangular aperture in one dimension illuminated by a monocromatic 
wave propagating in the z direction. The field distribution immediately after the aperture is given 
by 

(p(x,0) = l -> \x\<l/2, 

(p(x,0)=0 -> \x\>l/2, (27) 

with I being the aperture width. 
Considering that ip(x, 0) is real we can write 

4 + i)<4-i)= h (5 + I-m)h(1-£-|4 <28> 

We then apply Eq. (28) to the WDF definition Eq. (20) to find 

W{x,p x ) = -> \x\>l/2, 
W(x,p x ) = !1 -> 0<x<//2, 

W{x,p x ) = >± -Z/2<x<0, (29) 

After propagation we obtain the following integral field distribution 

\y( xz )\ 2 = ^l{f nx/z 2sm[kp x (l-2zp x /n-2x)} 

k 2 \Jn(2x-l)/(2z) kp x 

+ r^+O/W 2sm[kp x (l + 2zp x /n + 2x)] 1 

Jnx/z kp x \ 
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Fig. 1. Fresnel diffraction pattern for a one-dimensional aperture of width 0.1 mm with k = 10 7 . 

Fig. 1 shows a typical diffraction pattern obtained by numerical integration of Eq. (30). 

For wide angles paraxial approximation no longer applies and the appropriate WDF transformation 
is now given by 



W(q,p,z) = W(q 



n " ~ ipi 

W{q,p,z) = 



|P| < n, 
otherwise. 



(3D 



Eq. (31) shows that only the moments such that |p| < n can be propagated. 8 In fact, if |p|/n = 
sin a, with a the angle the ray makes with the z axis, It is obvious that the higher moments would 
correspond to values of | sin a \ > 1; these moments don't propagate and originate evanescent 
waves instead, Fig. 2. The net effect on near-field diffraction is that the high-frequency detail near 
the aperture is quickly reduced. 

The field intensity can now be evaluated by the expression 

k 2 \ip(x,z)\ 2 rpo 1 r r / i— 

- = J — sin \kp x \l-2\x- zp x / ^n 2 -pij^jdp x 



4vr s 



'pi kp x 

rV2 1 

+ / - 7 — sin \kp x 
Jpo kp x 



1 + 2 (x - ZPx/y/ri- 



-vl 



dp x , 



with 



Po 
Pi 



nx 



n{2x - I) 



{2x - I) 2 + Az 2 ' 



(32) 
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Fig. 3. Near-field diffraction pattern when the aperture width is exactly one wavelength; {k = 
10 6 m _1 , both scales in /jm). 

Fig. 3 shows the near-field diffraction pattern when the aperture is exactly one wavelength wide. 
The situation is such that all the high frequencies appear at values of \p x \ > 1 and are evanescent, 
resulting in a field pattern with one small minimum immediately after the aperture, after which the 
beam takes a quasi-gaussian shape, without further minima. The width of the sharp peak just after 
the aperture is considerably smaller than one wavelength determining a super-resolution on this 
region. 
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4 Special relativity 



Special relativity deals with a 4-dimensional space-time endowed with a pseudo-Euclidean metric 
which can be written as 

d s 2 = dx 2 + dy 2 + dz 2 - dt 2 , (34) 

where the space-time is referred to the coordinates (t, x, y, z). Here the units were chosen so as to 
make c = 1, c being the speed of light. 

For a more adequate optical interpretation one can use coordinates (x, y, z, r) with r the proper 



time 



.9 



T=t v-?. 



2\ V2 



t 
1 



Here v 2 = Ivl 2 , with v 



dx dy dz^ 
dt 1 dV dt 



(35) 



(36) 



the usual 3-velocity vector. Fig. 4 shows the coordinates of an event S when just two spatial coor- 
dinates are used together with coordinate r. 



Fig. 4. The relativistic frame in two dimensions. The r coordinate is the proper time while the 
distance to the origin is the time measured in the frame at rest. 

Considering Eqs. (35, 36) and the new coordinates the metric defined by Eq. (34) becomes 



ds' 



v 2 dt 2 - dt 2 = dr 2 



(37) 



The trajectory of an event in 4-space is known as its world line. In Fig. 5 we represent the world 
line of an event 8 with coordinates (x, r). At each position the derivative of the world line with 
respect to r is 

dx dt rrins 
x = — — — = ^v x . (38) 
dt dr 

More generally we can write 

7/ = TV, (39) 
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X 



Fig. 5. The curved line represents the word line of event £. The speed is represented by the vector 
v, which is the side of a rectangular triangle whose hypotenuse has a magnitude c. 



with 



v = {x,y, z 



(40) 



It follows from Eq. (40) that, at each point, the components of v are the direction cosines of the 
tangent vector to the world line through that point. 

We must now state a basic principle equivalent to Fermat's principle for ray propagation in optics 
given in Eq. (1). Taking into consideration the relativistic Lagrangian we can state the following 
variational principle: 3 ' 4 



where m is the rest mass of the particle and V is the local potential energy. It can be shown that 
TTry represents the kinetic energy of the particle. Using Eq. (37) a straightforward calculation shows 
that ds = dr so that we can also write 



The new principle means that the particle will choose a trajectory moving away from high potential 
energies much as the optical ray moves away from low refractive index. From this principle we can 
derive all the equations of relativistic mechanics, in a similar manner as from Fermat's principle it 
is possible to derive light propagation. 

Comparing this equation with Eq. (9), we can define the Lagrangian of the mechanical system as 



We now have a 4-dimensional space, while in the Hamiltonian formulation of optics we used 3 
dimensions. The following list shows the relationship between the two systems, where we refer 
first to the relativistic system and then to the optical system. 




(41) 




(42) 



L = (7717 — V) . 



(43) 



x <-> x, 



y <-> y, 



z <-> no equivalent, 



T <-» Z. 



In the mapping from optical propagation to special relativity the optical axis becomes the proper 
time axis and the ray direction cosines correspond to the components of the speed vector v. The 
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refractive index has no direct homologous; it will be shown that in special relativity we must 
consider a non-homogeneous medium with different refractive indices in the spatial and proper 
time directions. We can derive the conjugate momentum from the Lagrangian using the standard 
procedure: 

p = I I = mv. (44) 

Comparing with Eq. (13) it is clear that m is the analogous of the position dependent refractive 
index. 

The system Hamiltonian can be calculated 

H = p-v-L 

mv 2 

— — Lj 



7 
m 

= +V. (45) 

7 

The canonical equations follow directly from Eq. (15) 

It = TV, 

^ = -gradtf, (46) 

GST 

where the gradient is taken over the spatial coordinates, as usual. 

The first of the equations above is the same as Eq. (40), while the second one can be developed as 

d ( m_\ f m \ 

grad V . (47) 



dr \ t / \ T 

Considering that from the quantities inside the gradient only the potential energy should be a 
function of the spatial coordinates, we can simplify the second member: 

-gradl/ (48) 




= -gradV. (49) 

Eq. (49) is formally equivalent to the last two Eqs. (19), confirming that the conjugate momentum 
components are proportional to world line's direction cosines in 4-space. The total refractive index 
analogue can now be found to be m — V/j. We can check the validity of Eq. (49) by replacing the 
r derivative by a derivative with respect to t. 

— (rwyv) = — gradV, (50) 

(JjL 

where is the relativistic mass and the product rwyv is the relativistic momentum. 

If the mass is allowed to be coordinate dependent, as a mass distribution through the Universe, 
the passage between Eqs. (47) and (48) is illegitimate and we are led to equations similar to Eq. 
(19). The consideration of a coordinate dependent mass allows the prediction of worm tubes, the 
analogues of optical waveguides, and black holes, which would have optical analogues in high 
refractive index glass beads with gradual transition to vacuum refractive index. 
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5 De Broglie's wavelength 

The formal equivalence between light propagation and special relativity in the x, y, z, r frame 
suggests that de Broglie's wavelength may be the formal equivalent of light wavelength. We would 
like to associate an event's world line to a light ray and similarly we want to say that, in the absence 
of a potential, the event's world line is the normal to the wavefront at every point in 4-space. We 
must start from a basic principle, stating that each particle has an intrinsic frequency related to it's 
mass in a similar way as the wavelength of light is related to the refractive index; we state this 
principle by the equation 

(51) 



m 



where h is Planck's constant. If we remember that everything is normalized to the speed of light 
by c = 1, Eq. (51) is equivalent to a photon's energy equation 



E = hf. 



(52) 



So we have extended an existing principle to state that any particle has an intrinsic frequency that 
is the result of dividing it's equivalent energy E = mc 2 by Planck's constant. 

In the 4-space x, y, z, r frame a material particle travels in a trajectory with direction cosines 
given by the components of v and consistently a photon travels along a spatial direction with zero 
component in the r direction. The intrinsic frequency defined by Eq. (51) originates a wavelength 
along the world line given by 

c _ h 

A w — — — — , (53) 
/ mc 

where we have temporarily removed the c normalization for clarity reasons. 



T A 




Fig. 6. The moving particle has a world wavelength X w = h/(mc) and a spatial wavelength A = 

h/(mv). 

As shown in Fig. 6, when projected onto 3-space \ w defines a spatial wavelength such that 



A 



A,„c 



h 
V 



(54) 



The previous equation defines a spatial wavelength which is the same as was originally proposed 
by de Broglie in the non-relativistic limit. When the speed of light is approached Eq. (54) will 
produce a wavelength A — > \ w while de Broglie's predictions use the relativistic momentum and 
so A — > oo when v — > c. 
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6 Wave propagation and Schrodinger equation 

The arguments in the previous paragraphs lead us to state the general principle that a particle has 
an associated frequency given by Eq. (51) and travels on a world line through 4-space with the 
speed of light. In a generalization to wave packets, we will respect the formal similarity with light 
propagation and state that all waves travel in 4-space with the speed of light. A particle about 
which we know mass and speed but know nothing about its position will be represented by a 
monocromatic wave and a moving particle in general will be represented by a wave packet. 

According to general optical practice we will say that the field must verify the wave equation 

(v 2 - = 0, (55) 

where V is a point in 4-space and V 2 is an extended laplacian operator 



^ dx 2 dy 2 dz 2 ^ dr 2 ^ ^ 

In Eq. (55) we have returned to the c = 1 normalization in order to treat all the coordinates on an 
equal footing. Due to the special 4-space metric we will assume that (p(V, t) is of the form 

<p(V,t) = §(V)e i2nft , (57) 

with / given by Eq. (51). Notice that we used a plus sign in the exponent instead of the minus sign 
used in optical propagation; this is due to the special 4-space metric. 

Not surprisingly we will find that, in the absence of a potential, Eq. (55) can be written in the form 
of Helmoltz equation 

(V 2 + k 2 ) $(V) = 0, (58) 

with 

* = (59) 

If we take into consideration Eq. (35), the laplacian becomes 

„ d 2 d 2 d 2 2 d 2 



d 2 d 2 d 2 & 



dx 2 dy 2 dz 2 dt 2 
= + (60) 

where V 2 represents the usual laplacian operator in 3-space. 
In order to derive Schrodinger equation we re- write Eq. (55) 

VV(P,t)+zA3^^ = 0, (61) 

and using Eq. (60) 

V^tJ + ^f + l)^^. (62) 
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In a non-relativistic situation 7 
Schrodinger equation 10 



where h = h/(2-n). 

Eq. (63) retains the symbol V representing a point in 4-space. It must be noted, though, that in a 
non-relativistic situation r — > t and we can say that (p(V, t) — > £) with P having the 3-space 
coordinates of V. 

In the presence of a potential we have to consider that momentum is no longer preserved, as shown 
by Eq. (49). This can be taken into account when we evaluate the laplacian in Eq. (62) by the 
inclusion of an extra term —V(p(V,t). We will the end up with the Schrodinger equation in a 
potential. 



— > 1. Considering Eq. (53) we can write Eq. (62) in the form of 



7 Heisenberg's uncertainty principle 

For the pair of associated variables x and p x , Heisenberg's uncertainty principle states that there is 
an uncertainty governed by the relation 

AxA Px > ~. (64) 

The interpretation of the uncertainty principle is that the best we know a particle's position, the 
least we know about its momentum and vice-versa; the product of the position and momentum 
distribution widths is a small number multiplied by Planck's constant. An application of the uncer- 
tainty relationship is usually found in the diffraction of a particle by an aperture. 

If we assume that the localization of a particle with momentum p x can be done with an aperture of 
width Ax, we can use Fraunhofer diffraction theory to say that further along the way the probability 
of finding the particle with any particular value of it's momentum will be given by the square of the 
aperture Fourier transform, considering de Broglie's relationship for the translation of momentum 
into wavelength. 

A rectangular aperture of width Ax has a Fourier transform given by 

A(f x ) = Axsinc(Ax/ a; ), (65) 
where sine is the usual sin(x)/x function. 

Considering de Broglie's relationship given by Eq. (54), making f x = 1/A and the fact that the 
Fourier transform must be squared we can write 

Pip.) = Ax 2 sine 2 (^p) . (66) 

The second member on the previous equation has its first minimum for Axp x /h = n and so we 
can say that the spread in momentum is governed by AxAp x = irh and Eq. (64) is verified. 
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If we accept that wave packets propagate in 4-space at the speed of light and that the momentum 
is given by Eq. (44), there is a upper limit to the modulus of the momentum |p| < mc. In the 
propagation of light rays we found a similar limitation as |p| < n with the results in light diffraction 
exemplified by Eq. (32) and Fig. 3. It is expected that the same effects will be present in particle 
diffraction and in fact Figs. 2 and 3 could also represent the diffraction of a stationary particle by 
an aperture with width equal to A„,. The strong peak about half one wavelength in front of the 
aperture shows that the particle is localized in a region considerably smaller than its wavelength 
and, above all, shows the inexistence of higher order peaks. 

8 Conclusions 

Special relativity was shown to be formally equivalent to light propagation, provided the time axis 
is replaced by the proper time. In this coordinate set all particles follow w world line at the speed 
of light and can be assumed to have an intrinsic frequency given by mc 2 /h. Quantum mechanics is 
then a projection of 4-space wave propagation into 3-space. Important conclusions were possible 
through the analogy with light propagation and diffraction. 

It was possible to derive Schrodinger equation and it was shown that Heisenberg's uncertainty 
principle may be violated in special cases in the very close range, similarly to what had already 
been shown to happen in light diffraction. 1 

Future work will probably allow the derivation of relativistic quantum mechanics conclusions, 
through the use of the Wigner Distribution Function for the prediction of wave packet propagation 
in 4-space. 
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